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We elaborate on a general method to find complex envelope solitons in a cigar shaped Bose-Einstein 
condensate in a trap. The procedure incorporates time dependent scattering length, oscillator fre- 
quency and loss/gain. A variety of time dependencies of the above parameters, akin to the ones 
occurring in the experiments can be tackled. 



Coherent atom optics is a subject of much current interest for its relevance to both basic 
physics and technology. For this purpose, lower dimensional condensates e.g., cigar shaped 
Bose-Einstein condensates (BECs) have been the subject of active study in recent years 
PQ. Intense investigations about the behavior of the condensate in the presence of time 
varying control parameters like, scattering length, oscillator frequency etc., are being carried 
out for the purpose of optimal control. In BEC, solitary waves have been experimentally 
observed, both in repulsive and attractive domains [21 E]- These have been understood from 
the non-linear Gross-Pitaevskii (GP) equation describing weakly coupled BEC, 

m ^&r = {-^ + UMr,t)\* + v}*(T,t), (1) 

where U = 4irh 2 a s (t)/m, a s = the s-wave scattering length, and V is the external potential. 
The scattering length, which is the coefficient of nonlinearity in GP equation, can be adjusted 
both in sign and magnitude through Feshbach resonance. Recently a family of exact solutions 
of quasi-one dimensional GP equation with time varying scattering length, loss/gain, in the 
presence of oscillator has been obtained [I]. These are chirped solitons with real envelope. 
In recent years, solitons with complex envelope is generating interest in different areas of 
physics. These are known as Bloch solitons in condensed matter physics. In optical fibers 
one also finds these solitons. In BEC, complex solitons have been investigated (5j El [Zl [HI IS] - 
Very recently complex envelope solitons in optical lattices has been reported |lUj . 

In the present study, we explicate a method to obtain general complex envelope solitons 
in cigar shaped BEC, in a trap incorporating time dependent scattering length, oscillator 
frequency, and loss/gain. We have constructed the two hydrodynamic equations in this 
scenario and have found their solutions. Our analysis starts with quasi-one dimensional non- 
linear Schrodinger equation (NLSE) derived from three dimensional GP equation with an 
additional loss/gain term g(t) and time dependent chemical potential f(t). In dimensionless 
units this is given by [I], 

idtf = -^ + iM(t)^ + 7 (t)|^| 2 ^ + ^^ + ^W- (2) 
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We assume the following ansatz solution 

i/>(z,t) = J~Mt)F[A(t)z - Z(t)]exp 



i{x(z,t) + <f>{z,t)} 



G(t) 



(3) 



where \jA(t)F[A(t)z — l(t)] exp[ix(z, t)] carries the signature of complex envelope. Further 
G(t) and l(t) have been defined as, 

G(t) = [ g(t')dt' , l(t)= [ v(t')dt'. (4) 
Jo Jo 

The phase is chirped as, 

1 lr* 

(f)(z,t) = ait) + b(t)z - -c(t)z 2 , where, a(t) = a - - A 2 (t')dt'. 

2 2 jo 

Substitution of the above ansatz in Eq.(T5]) allows one to separate this equation in imaginary 
and real parts. The knowledge of the earlier conditions, leads one to write the imaginary 
equation as, 



-uF' = -F'x' - ^Fx", and 



X 



u — 



2Co 
F 2 



(5) 



Here prime denotes derivative with respect to T, where T = ^4(t)[z — l{t)]. We have set 
it+ci-A = u ^ A t = j± c anc j q o j s an integration constant. Redefining Co = ^if-, F = y/a and 



cr , Eq.(jS]) leads to, 



X' = u(l | . 



a J 



(6) 



We have imposed the boundary condition, x' ~^ for a — > gq. The local phase velocity 
of the solitary pulse is described by x' an d local density is described by g, g represents 
equilibrium density. 

From the real equation one can extract a Riccati type equation, c t — c 2 (t) = M(t) which 
yields the solution of c(t). The real equation also yields, 



F" -uF- 2kF 3 



X ' 2 F - 2u X 'F, 



(7) 



where k — Applying the following consistency conditions, one obtains 



7 (t) = lQ e~ G A{t)/A , b(t) = A(t), 



A(t) 



A exp 



c(t')dt' 



, v(t) = A 2 ^ 



and, 



A 



F" -eF- 2kF a = 

F 6 



where e = fi — u 2 and A = u 2 Fq. Integration of Eq.Q leads to the convenient form, 

2 / 



[KG — U 



(G - G ) 



G 



(9) 



The solution for Eq.flH]) is, 

a{z, t) = ctq < 1 — cos 2 9 sech 2 



A(z - I) cos 9 



c 



(10) 



These are the dark and grey solitons in the distributed scenario. 
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Figure 1: (a) Density profile of the Condensate and (b) Chirp Controlled Soliton Motion 

In conclusion, for quasi-one dimensional GP equation in an oscillator potential with time 
dependent coupling and loss/gain, a general complex dark soliton profile is obtained. In the 
presence of expulsive oscillator the motion of the condensate is controlled by the sign of the 
chirp, i.e, whether the soliton is accelerating or decelerating. 
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